SOS model of overlayer induced faceting 
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A solid-on-solid model is proposed to describe faceting of bcc(lll) metal surface induced by a 
metal overlayer. It is shown that the first order phase transition occurs between faceted {211} or 
{110} and disordered phases. The ordered phases consist of large 3-sided pyramids with {211} facets 
or {110} facets. It is shown that the high-temperature disordered phase has not planar bcc(lll) 
structure but faceted disordered structure. Hysteresis effects were observed when the system was 
warmed above the transition temperature and then cooled down. Temperature dependence of LEED 
patterns for faceted and disordered phase are calculated in kinematic approximation. 

PACS numbers: 68.35.Rh, 68.43.De, 64.60.Cn 
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I. INTRODUCTION 

Recent experiment»i*2i2i2i£i£i£ have demonstrated that 
surfaces such as W(lll) and Mo(lll) covered by a single 
physical monolayer of certain metal, e.g., Pd, Pt, undergo 
massive reconstruction from planar morphology to micro- 
scopically faceted surface after annealing to T > 700A'. 
The reconstructed surface consists of 3-sided pyramids 
with mainly {211} facets, and pyramid dimensions range 
from ~ 1 to 100 nm. It has been shown that facets are 
composed of substrate atoms, and the monolayer of ad- 
sorbate remains on the outermost surface layer during 
the faceting transformation. Another type of massive re- 
construction has been very recently found in STM and 
LEED experiment for Pd on Ta(lll) system^. The recon- 
structed surface consists of {011} facets which form large 
triangular pyramids. The third type of reconstructed 
structure, coexistence of small {011} facets with large 
{211} facets, has been observed 2 in Pd on W(lll) for 
prolonged annealing time (for short annealing time only 
{211} pyramids occurred). An important step in under- 
standing the thermal stability of reconstructed surfaces 
was LEED experiment performed in high temperatures 
for Pd on Mo(lll) by Song et ali They demonstrated 
existence of reversible phase transition between faceted 
{211} and planar phases and they found that this tran- 
sition has a large hysteresis, i.e., the transition tempera- 
ture in cooling cycle (T ss 830) is lower than that in the 
heating cycle (T w 870). 

In theoretical studies of overlayer- induced faceting 2 **^!! 
the first principle method have been used to calculate 
the surface formation energy of fee metals (Pd, Pt, Au, 
Ag, Cu) adsorbed on Mo and W. Results of these calcu- 
lation performed for pseudomorphic adsorbate overlayer 
on (111), (211), and (011) flat surfaces show that (111) 
surface becomes unstable at the coverage of one physical 
monolayer, where the energy of (211) orientation has the 
lowest value. However, such energy calculation (at T=0) 
are not sufficient to explain why some of these metals 
do not induce faceting, e.g., Ag/Mo(lll). On the other 
hand, the first principle calculation confirmed that criti- 
cal coverage to induce faceting is approximately equal to 



one physical monolayer (PML) and that the growth mode 
is Stransky-Krastanov as the surface energy increases at 
coverage higher than 1 PML. 

Theoretical studies of surface reconstruction and sur- 
face phase transitions in bimctalic system are com- 
plicated problems mainly due to occurrence of long 
range many-body interactions. One of approaches to 
such problems is application of simple solid-on-solid 



model; 
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in which the crystal is represented 



by two-dimensional array of columns. They were em- 
ployed to study roughening transition 1112 , missing- 
row reconstruction^, growth of the surface^, surface 
diffusioni 2 * 1 ^ adsorption 14 , phase transitions in two com- 
ponent crystal^, etc. 

In this paper we introduce a simple SOS model for 
bcc(lll) surface covered by one physical monolayer of 
adsorbed atoms to describe surface reconstruction and 
phase transitions. Using Monte Carlo simulation we 
study change of the surface structure during the heat- 
ing and cooling processes, phase transitions to a faceted 
phases, phase diagram, LEED patterns and temperature 
dependence of diffracted intensity calculated in kinematic 
approximation. 



II. THE SOS MODEL 

To study faceting in bimetalic system one needs to 
know interaction potentials. There are many-body po- 
tential for Mo-Mo and W-W interactions derived by Fin- 
nis and Sinclai»i£ but to our knowledge there is only Ni- 
Mo many-body potential^ for interaction between (Mo, 
W) and fee metals. However in Ni/Mo(lll) system the 
faceting does not occur—. Therefore we are going to study 
overlayer induced faceting using a simple solid-on-solid 
(SOS) model which describe a surface formation energy 
in agreement with result of the first principle calculation. 

In order to construct a model of overlayer-induced 
faceting on bcc(lll) surface (BCCSOS) we are taking 
into account some experimental evidences. There is a 
critical coverage, approximately equal to 1PML, to in- 
duce faceting on W(lll) and Mo(lll)^. When the cover- 
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FIG. 1: Schematic view of the bcc(lll) surface. Sites of 
the sublattices / = 0,1,2 are denoted by squares, triangles, 
crosses, respectively. The Z axis is normal to the plane and 
a x = ay/2/2, a y = aVE/6. 



age exceeds 1PML, extra adsorbate atoms form 3D clus- 
ters on the top of faceted surfaces. Hence to simplify 
problem it will be assumed that coverage is constant and 
equal to 1PML. Such coverage is equivalent to 3, 2, and 
1 geometric monolayers of adsorbate on bcc (111), (211), 
and (Oil) surfaces, respectively. A second assumption is 
that substrate atoms (B) can take position in discrete bcc 
lattice. It has been shown£, in the case of Pd/W(lll), 
that both {011} facets and {112} facets of W are cov- 
ered by pseudomorphic monolayer of Pd. Therefore we 
assume that adsorbate atoms (A) are also located in 
positions of the bcc structure. There are three types 
of interaction energies in this SOS model: adsorbate- 



adsorbate interaction e 



AA 



adsorbate-substrate interac- 



tion ef B , and substrate-substrate interaction e BB , with 
the range of interaction up to forth neighbors (i =1, 2, 3, 
4). 

It is convenient to chose the coordinate system with the 
Z axis parallel to [111] direction and X, Y axes along 
[101], [121] directions in the (111) plane. Atoms along 
the closed packed direction ( parallel to the z axis) form 
columns which positions in the (x, y) plane are described 
by three triangular sublattices I = 0, 1, 2 shown in Fig.^ 
Column height %i at site i in the Ith sublattice is defined 
as Zi = ft, J av / 3/6 where a is bcc lattice constant, hi = 
3rii + 1 and is the number of atoms in ith column. The 
assumption on constant coverage equal to 1PML means 
that there is exactly one A atom in each column placed 
on its top. As we are going to study reconstruction only, 
not desorption or crystal melting, the typical restriction 
of SOS models on columns heights will be assumed. The 
hi difference between nearest neighbor sites are forced to 
be ±1,±2. 

Let us define the Hamiltonian as surface formation en- 
ergy of A material on B (111) surface 



where Ea/b(Na, N b ) denotes energy of Na atoms A on 
bcc(lll) surface consisting of N B atoms B, E A ulk (NA), 
E h g lk (NB) are energies of A, B atoms in their bulk en- 
vironments. Now taking into account above assumptions 
we can express Hamiltonian Eq. Q in terms of column 
heights hi in the following form 



H = \ E E [Ji5(\lH - h n | - 1) + K,S (\hi - h n | - 2)] 

I [ 31 

+ £ [2J 2 S(\hi - h h |) + (2 J 2 + K 2 ) 6 (\hi - h h | - 3)] 

32 
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and sums over ji, j 2 , and js, denotes summing over first, 
second, and third neighbors of a column at site i. In 
what follows, we will treat J\, J 2l K%, and K 2 as model 
parameters. It is worth noting that this BCCSOS Hamil- 
tonian of overlayer-induced faceting at constant coverage 
can also be used to study a clean bcc(lll) surface by 
setting ef B = e 



AA 



r BB 



A. Energies of (111), (211), and (110) faces at 
T=0 



First, we will check the stability of ideal surfaces (111), 
(211), and (110) covered with 3, 2, 1 geometric monolay- 
ers of A, similarly as it was performed in first princi- 
ples calculations^. This allows us to estimate values of 
model parameters for reconstruction, from (111) surface 
to {112} or to {011} faceted surface, under assumption 
that edge energies are neglected. 

Using the Hamiltonian Eq. we get the following 
expressions for surface energy per site for ideal face of 
orientation (hkl) 



H = E a/ b(Na,N b ) 



Ea 



E b B ulk (N B ), (1) 



Em = E, 
Eno = Em 4 
E 2 u = Em 4 
where 

E r = J + if 1 



7J 2 , 
f 2J 2 
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Thus, the stability conditions are: 

K 2 > 0, J 2 + K 2 > for (111) surface, 

J 2 + K 2 < 0, 2J 2 +K 2 <0 for (110) surface, 

K 2 < 0, 2J 2 + K 2 > for (211) surface. 

It easy to see that transformation from (111) surface to 
{211} facets will be possible for negative values of K 2 and 
appropriate positive J 2 . To satisfy these conditions, the 
range of interactions e" 13 should be not smaller than — 
a\/TT/2. On the other hand, it is easy to ensure stability 
of bcc(lll) surface at zero coverage ( ej = 0) by choosing 
e B B < and e BB < 0. Thus, there is general possibility 
to chose a set of interactions {e^ B , e BB , €3, £4} in such a 
way that at zero coverage the bcc(lll) surface is stable 
whereas for coverage equal to one physical monolayer the 
bcc(211) surface become stable. 

It is worth noting that surface energies per site Ehki 
are related to surface energies per surface atom OhM (used 
in first principles calculation^) in the following way 
am = 3-Em, (72ii = j Clio = ^HOj as the number 

of surface atoms is N/3, N/2, and N, respectively. More- 
over, above stability conditions correspond to faceting 
conditions^ expressed in terms of scaled <Jhk.li e.g., the 
conditions E 2 \\ < Em, E 2 \\ < £110 are equivalent 
to 3(72ii/2 < (Tin, 3(72ii/2 < 3(7iio for transformation 
(111) -{112}. 

Independence of stability condition on parameters Jo, 
Ji, and K\ comes from assumptions on constant num- 
ber of adsorbate atoms A and the restriction on column 
heights. Therefore we neglect these parameters in further 
analysis what corresponds to shifting of the surface en- 
ergy by — E r . Moreover, by choosing J 2 as the unit of en- 
ergy we will work with dimcnsionlcss quantities: energy 
H = H I ' J 2 , temperature T = ksT / J 2 , and parameter 
K = K 2 /J 2 (in what follows, the tilde will be omitted). 

III. SIMULATION OF FACETING 

We will discus here results of the annealing process in- 
vestigated by Monte Carlo (MC) simulation in canonical 
ensemble. Let us start with some details concerning the 
simulation method. Configurations of BCCSOS model 
are generated via the classic Metropolis algorithm. A 
new configuration in Markov chain is generated from the 
previous configuration by moving one B atom from a site 
! to a site j, what is equivalent to the following change 
of heights: (hi,hj) — > (hi — 3,/ij + 3). It is important 
to note that such change of heights is possible only if 
hi is the local maximum ( with respect to the nearest 
neighbors) and hj is the local minimum because of the 
restriction on column heights. Let us remember one of 
the model assumption that adsorbate atoms always stay 
on the column tops as we study reconstruction of B sur- 
face covered by the physical monolayer of A material. 
MC simulation were carried out on rectangular lattice of 
linear size L x ,L y along x, y direction, respectively with 
periodic boundary conditions (PBC). To check the role of 
PBC additional simulations were performed with bound- 
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FIG. 2: Temperature dependence of the surface energy dur- 
ing the warming up cycles (triangles) and the cooling cycles 
(circles). 



ary atoms fixed at positions of the flat bcc(lll) surface. 
The results are consistent with those obtained for calcu- 
lation with PBC. 

First, we investigate properties of the system with cou- 
pling constant K = —1.25 during warming up process. 
At each temperature Tj = + AT system spends the 
same annealing time r measured in Monte Carlo steps 
per site. The ideal bcc(lll) surface is used as a start- 
ing configuration at To. We also investigate properties 
of the system during cooling down the sample. Results 
of simulations for the following parameters: Lx = 192, 
Ly = 336, AT = 0.03 and t = 3x 10 5 are presented in 
Fig. HJ- Fig. El Let us notice that these sizes of the lattice 
correspond to the area of 426A x 43lA on the bcc(lll) 
surface with the lattice constant a = 3.14A. 

The surface formation energy decreases as tempera- 
ture is elevated (see Fig. what indicates reconstruc- 
tion of the surface. At T w 0.4 very small 3-sided {211} 
pyramids are formed on the surface and further warm- 
ing up causes growth of pyramids sizes (see Fig- EJ - By 
the {211} pyramid we mean the pyramid built of (211), 
(121), and (112) facets. A big jump of the surface energy 
is observed at T — 1.26 where the phase transition to dis- 
ordered phase occurs. On the other hand, discontinuous 
change of the surface energy is observed at T = 1.20 when 
the system is cooling down. The hysteresis loop seen in 
Fig. |21 indicates that the phase transition is of first order. 
The surface is undercooled or overheated because system 
rest in metastable state separated from the stable state 
by the free energy barrier. The hysteresis effects, depen- 
dence of temperature of the phase transition on cooling 
and warming, has been observed in LEED experiment'' 
for Pd/Mo(lll). When system is cooling down we ob- 
serve formation of large pyramids with defected faces and 
edges just below the temperature of the phase transition 
(see Fig. 0J. Facets of similar shapes and sizes were ob- 
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FIG. 3: Snapshot of the surface at T=0.72 on warming. 




FIG. 4: Snapshot of the surface at T=1.17 on cooling. 



served at this temperature on warming. During further 
cooling defects disappear and pyramids take nearly ideal 
shapes at T w 1. Therefore the dependence of surface 
energy on cooling is very week in the ordered faceted 
phase. 

It is very interesting that in the disordered phase the 
surface is not planar as suggested in experimental works 7 . 
The structure of disordered phase contains randomly dis- 
tributed small facets mainly of {112} orientations what is 
shown in Fig. [5] where the cross section of the disordered 
phase is compared to cross sections of faceted {211} and 
ideal bcc(lll) surfaces. Moreover, we calculated the av- 
erage number of atoms on a facet of orientation (211) 
which in the disordered phase is nearly 10. So we will 
call this phase as disordered faceted phase (DFP). In very 
narrow temperature range, we observed a coexistence of 
a faceted {112} and DFP, e.g., in simulation at T = 1.23 
denoted by the triangle in Fig. |2 A surface structure 
in this case looks like surface of ordered phase where 
large pyramids are removed and on their place DFP is 
present. On the other hand, there are 3-sided pyramidal 
holes built of large facets of {112} orientation. Some- 
times, DFP is located also at the bottom of these holes. 
It is worth noting that coexistence of both faceted and 
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FIG. 5: Plot of column height along the y direction for the 
surface: (a) below the phase transition, T=1.17, (b) above 
the phase transition, T=1.50, (c) ideal bc(lll) - initial con- 
figuration. 



disordered phases has been observed^ in LEED experi- 
ment in Pd/Mo(lll) system. However, it was suggested 
that coexistence is due to the inhomogeneity of the sur- 
face. 

One of quantities measured in MC simulation is the 
square mean width of the surface, 



Sh 2 = 



(h)f 



(3) 



Behaviour of the Sh 2 during warming up and cooling 
down processes allows us to study dependence of mean 
vertical sizes of pyramids on temperature. In Fig. 
we present results of simulation for 3 different anneal- 
ing times r = 2 x 10 3 , 2 x 10 4 , 3 x 10 5 . In all cases 
the hysteresis loop is present near the phase transition 
and for shorter time a larger hysteresis is observed. It is 
easy to explain because it is more probable to overcome 
the free energy barrier in longer time. In the warming 
up process the size of pyramids does not depends on the 
annealing time r up to T rs 0.5. When temperature ap- 
proaches T w 1 then a rapid growth is observed and for 
longer annealing time r greater pyramids are formed on 
the surface. Pyramids reach maximal sizes just below the 
temperature of the phase transition where dependence of 
sizes on annealing time is the strongest. It worth noting 
that dependence of facets size on temperature qualita- 
tively agrees with the experimental results^. The value 
of Sh 2 rapidly decreases above the phase transition where 
very weak dependence on temperature is observed. In 
the cooling down process large pyramids are formed just 
below the phase transition temperature, and further de- 
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FIG. 6: Temperature dependence of the Sh 2 during the heat- 
ing and the cooling cycles for 3 different annealing times. 
Each label contains symbols for heating, cooling and the num- 
ber of MC steps 



creasing of temperature practically does not change the 
size of pyramids. 

We check dependence of surface growth on the initial 
temperature To and we find that the behaviour of Sh 2 
near the transition point is very similar to that in Fig. 
for process which starts at To = 0.3 or To = 1-0. Finally 
we analyze the influence of the temperature increment, 
AT, on the surface energy and the Sh 2 near the phase 
transition for r = 3 x 10 5 . We found that for AT = 
0.01 the hysteresis loop slightly increases, there are more 
points where the faceted and disordered phases coexist 
but it seems that sizes of pyramids do not depend on the 
temperature increment. 

We observe that the growth of the facets on warm- 
ing base on formation of a larger pyramid from several 
smaller pyramids. Probability of such rebuilding of the 
surface is very small at low temperature. This is reflected 
in the acceptance ratio /, a fraction of accepted config- 
urations in the Markov chain, which is very low at low 
temperatures, / « 0.002 below T = 0.3 and / = 0.01 
at T = 0.7. On the other hand, the acceptance ratio 
is of the same order for T = 1.0, / = 0.1, and close to 
the phase transition, / = 0.2, but we observe difference 
in facets sizes at these temperatures. Looking at snap- 
shots of the surface after annealing at temperature up 
to T s» 1.1 on warming, we found that pyramids have 
nearly ideal shapes (see for example Fig. |3J). We think 
that all above observation can be explained in the follow- 
ing way. The system is passing through metastable states 
- local minima of the free energy - as temperature is in- 
creased. There are barriers between local minima hence 
the growth of facets depends on the annealing time. 
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FIG. 7: Phase diagram: T versus K. 



IV. PHASE DIAGRAM 

The stability anal ysis of flat (111), (211), and (110) 
surfaces in Sec. Ill Al shows that at T = the (211) sur- 
face is stable for — 2 < K < whereas for K < —2 the 
(110) surface is stable. Whether these findings are gen- 
erally correct we can check constructing phase diagram 
in the (T, K) plane by using MC simulations. Most of 
simulations were performed at constant K varying tem- 
perature. To check the results some simulations were 
performed at constant temperature. To identify differ- 
ent phases we calculated masses of clusters of different 
facet orientations (111), {HO}, and {112}. A cluster of 
orientation (hkl) is defined in similar way as percolation 
cluster on the (hkl) plane. We also recorded equilibrium 
configurations at each investigated point of the (T, K) 
plane. We found that for K < the phase diagram ( see 
Fig. [7J contains three phases 

1. Faceted {112} phase consisting of 3-sided {112} 
pyramids. 

2. Faceted {110} phase consisting of 3-sided {110} 
pyramids. 

3. Disordered faceted phase. 

The surface of the faceted {110} phase is built of (HO), 
(101), and (011) facets which form 3-sided pyramids on 
the bcc(lll) surface. Such surface has been observed ex- 
perimentally in Pd on Ta(lll) system^. The phase tran- 
sition between the disordered phase and one of ordered 
is of first order as surface energy, the square mean width 
of the surface Sh 2 , and the structure factor (see Sec. 
|Vj change discontinuously at the transition point. The 
transition temperature between F{211} and disordered 
phase is bounded from above by maximum T c = 1.20 at 
K = —1.50. If we chose this point to estimate model 
parameters for Pd/Mo(lll) with T c = 850 K we get 
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J 2 = 0.06 eV, and K 2 = -0.09 eV. On the other hand 
the transition temperature from F{110} to the disordered 
phases is increasing function of \K | and it might reach 
large value. The present model does not include the 
desorption, thus we can expect that transition tempera- 
ture might be above the desorption temperature in some 
systems especially where the faceted F{110} phase oc- 
curs. In such cases the phase transition to the disordered 
faceted phase could not be observed but rather decon- 
struction induced by desorption would be expected. 

There is difficult to study phase transition between or- 
dered phases using MC simulation because for K w —2 
we observed coexistence of two types of facets {011} and 
{112}. The mixed phase does not contain separated pyra- 
mids of two types but mixture of facets, e.g., on the (211) 
facet a smaller (011) facet can occur. 



V. SIMULATION OF LEED PATTERNS 

Low energy electron diffraction (LEED) experiments 
play an important role in investigation of structures 
and phase transitions of faceted surfaces^. Very re- 
cently, Song et al, 7 where able to study LEED pat- 
terns of Pd/Mo(lll) system in high temperatures up 
to 1200K. They demonstrated the existence of reversible 
planar /faceted phase transition with transition temper- 
ature different for heating cycles (T rj 870A") and for 
cooling cycles (T w 830X). 

In this section we employee the BCCSOS model to 
analyze the temperature dependence of LEED patterns. 
The diffracted intensity is proportional, in the kinematic 
approximation, to the structure factor 



S(k) = 



^exp(ikRj) 



(4) 



where summing is over surface atoms, k = k/ — k^, and 
ki, k f is the wave vector of incident electron, scattered 
electron, respectively. Position of surface atoms is ex- 
pressed by column height h rj in site Tj = (Xj,yj) in the 
following way R,, = a(2>\/2xj, v6yj, v3h rj )/6. In a SOS 
model not all atoms from column tops are surface atoms. 
Therefore in calculation of diffracted intensity sometimes 
the shadowing factor is usedi^. In this paper we consider 
the contribution to S only from surface atoms. As a sur- 
face atom in any configuration we choose the atom from 
the top of column which has at least four lower nearest 
neighbors. In the case when number of higher columns is 
equal to 3 then atom is regarded as surface atom if is not 
surrounded by the higher neighbors. Such case occurs for 
atoms on {110} facets. We will calculate structure factor 
only for a set of k/, chosen in such a way to describe 
diffraction on faces (111) and {211}. For each of these 
faces we construct a pair of inverse lattice vectors gi, g 2 
parallel to the face, e.g., for (121) gi = 2tc/cl(\/2/2, 0, 0) 



(a) 




(b) 



FIG. 8: LEED patterns of (a) disordered phase where spots 
are represented by squares and also spots of ideal bcc(lll) 
surface (crosses) are shown, (b) faceted {211} phase where 
circles, diamonds, and triangles denote spots coming from 
diffraction on facets (121), (211), and (112), respectively. 



and g 2 = 271-/0(0,4^/9,2^3/9). Then, the wave vec- 
tors of scattered electrons, kf , for given face are chosen to 
satisfy condition of constructive interference (kj — kAi = 
Zgi + mg 2 , where l,m = 0, ±1, . . . and the symbol j| de- 
notes vector component parallel to the face. 

We calculated structure factors for the wave vector of 
the incident electron, k{ = 27r/a(0, 0, 1.43), normal to 
the bcc(lll) surface. Thus energy of incident electrons 
is equal to 31 eV for a = 3.14A, similarly as in LEED 
experiment for Pd/Mo(lll)t Figure|S|shows LEED pat- 
terns obtained in computer simulation for the model pa- 
rameter K = —1.25. In the ordered phase {211} we see 
besides 3 clovers (groups of 3 spots) observed experimen- 
tally, the forth clover in the center of the image, and six 
single spots more distant from the center. The pattern 
of disordered phase consists of three spots which posi- 
tions correspond to the centers of external clovers in the 
diffracted pattern of the faceted phase. Apparent lack 
of spots from bcc(lll) surface confirms that disordered 
phase has not the planar structure. 

Temperature dependence of the structure factors were 
calculated using MC simulation and changing tempera- 
ture in the way described in Sec. IIIII Most of simulations 
were performed on the lattice with Lx = 90, Ly = 132. 
The spots of disordered phase have the same diffracted 
intensities. Their temperature dependence (see Fig. [5]) on 
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FIG. 9: Temperature dependence of the structure factor for 
the wave vector corresponding to a spot of the disordered 
phase on warming (triangles) and on cooling (squares). The 
hysteresis loop for large lattice Lx = 192, Ly = 336, is shown 
in the inset. 



warming and on cooling shows the presence of hysteresis 
loop in behaviour of S(T) close to the phase transition 
what is in agreement with experimental results 7 . The 
hysteresis loop calculated for large lattice with Lx = 192, 
Ly = 336 is shown in the inset in Fig. El There are several 
points where the structure factor takes intermediate val- 
ues because at these temperatures coexistence of faceted 
and disordered phases has been observed. More coexis- 
tence points where found on warming than on cooling. 
We observe also discontinuity of the structure factor at 
the transition temperature. Thus hysteresis and discon- 
tinuity of the structure factors are another argument for 
the first-order phase transition between faceted and dis- 
ordered phases. During the process of warming up the 
system the S decreases with temperature and becomes 
very small when on the surface large enough {211} facets 
are formed. 

The diffracted intensities of a spot in the ordered phase 
was calculated for k / fulfilling the condition of construc- 
tive interference and also for vectors from a very small 
surroundings of k/ . Then maximal value of S were cho- 
sen to represent the intensity of spot at this ky. Wc 
use such procedure because the surface does not consist 
of ideal flat {211} facets during warming up and cool- 
ing processes, therefore spots can change shapes, intensi- 
ties, and positions. The diffracted intensities of spots of 
faceted phase are different because areas of facets of dif- 
ferent orientation are not equal. Moreover, their depen- 
dence on temperature during warming the system shows 
irregular behaviour (see Fig. ITU|) because of facets grow 
i.e., from several small facets a bigger one arises. We 
observe rapid growth of S at temperature T rs 0.9 when 
large pyramids are formed on the surface. During the 
cooling cycles, the S reaches nearly maximal value just 
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FIG. 10: Temperature dependence of the structure factor for 
the wave vector corresponding to a leaf of the right clover on 
warming (triangles) and on cooling (squares). 



below transition temperature what confirms formation 
of large {211} facets. The maximal values of S in the 
faceted phase are nearly ten times smaller than in the 
disordered phase because in former case only ~ 1/3 sur- 
face atoms contribute to S, i.e., surface atoms placed on 
facets of the same orientation. 



VI. DISCUSSION 

We have presented here the BCCSOS model to study 
overlayer induced faceting. Although it might seem that 
model is very simple it gives many results in agreement 
with experiments. It has been shown that bcc(lll) sur- 
face covered by a physical monolayer undergoes recon- 
struction upon annealing. We obtained two types of re- 
constructed surfaces: faceted {211} covered by 3-sided 
pyramids as observed in Pd/Mo(lll)£ and faceted {011} 
surface found in Pd/Ta(lll)&. The sizes of facets depend 
on annealing temperature and annealing time. At high 
temperature the phase transition to disordered phase oc- 
curs. We have shown that LEED diffraction patterns 
can be calculated, in kinematic approximation, in warm- 
ing up and cooling down processes for different electron 
energies. In dependence of diffracted intensities on tem- 
peratures the hysteresis effect was observed close to tem- 
perature of the phase transition. 

On the other hand, we can investigate some prop- 
erties not reported in experimental papers. First of 
all the phase transition can be studied in details using 
this model. It is shown that the phase transition is of 
first order because quantities such as surface energy, the 
structure factor change discontinuously as temperature 
reaches the transition point. Moreover the hysteresis ef- 
fects are observed when system is warming up and then 
cooling down. We can analyze the structure of the surface 
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in disordered phase as well as in the coexistence of the 
faceted and disorder phase close to the transition tem- 
perature. We found that the disordered phase consists of 
many small {211} facets randomly distributed. Hence it 
is not planar bcc(lll) surface. It would be interesting to 
check experimentally this prediction. It was suggested^ 
that coexistence of faceted and disordered phases is due 
to the inhomogeneity of the adsorbate coverage but we 
have demonstrated that the coexistence could also exist 
for homogeneous coverage. 

It seems that this model can be applied to study recon- 
struction of curved bcc surfaces where reconstructed sur- 
face has different structure than faceted bcc(lll) surface. 
In case of Pd deposited on a needle-shaped tungsten^, 
step like {211} microfacets has been observed. Using the 
present model we have not obtained the faceted phase 
with large {211} pyramids and small {011} pyramids as 



observed in Pd/W(lll)£ for prolonged annealing time. 
This might be due to lack of density fluctuation. So we 
are going to extend the model including dependence on 
coverage. The most general extension should base on re- 
placing of the interactions e" by many-body potentials. 
To construct such potentials, e.g., for Mo-Pd, W-Pd sys- 
tems, the results of the first-principle calculations could 
be usedi^. 
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